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A conventional approach to precision calculations of Higgs boson observables uses 
quark masses nic and ruf, as inputs. However, quark masses are single numbers 
that hide a variety of low-energy data from which they are extracted, and also 
hide the various sources of theoretical uncertainties and correlations with additional 
input parameters such as Higher-precision calculations, which are needed to 
give meaning to future measurements, require more direct engagement with the low- 
energy data in a global analysis. We present an initial calculation in this direction, 
which illustrates the procedure and reveals some of the theory uncertainties that 
challenge subpercent determinations of Higgs boson partial widths. 
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I. INTRODUCTION 

The discovery of the Higgs boson DEI marks the beginning of a new era for precision 
studies. Not only is unprecedented precision achieved in Standard Model (SM) calcula¬ 
tions nn with the knowledge of the Higgs boson mass miH], but experimental data on 
a large number of Higgs observables [U] allows us for the first time to scrutinize the Higgs 
sector of the SM m and beyond niHn]. Any discrepancy between precision data and SM 
predictions would be an indication of new physics. 

Though not explicitly stated in the context of precision Higgs analysis, an important role 
in this program is played by low-energy observables, such as moments of e^e“ annihilation 
cross section and moments of semileptonic B decay distributions. In fact, our knowledge of 
the charm and bottom quark masses tuq {Q = c,b), which are important inputs of precision 
Higgs calculations, largely comes from analyzing these low-energy data. This can be seen 
from the fact that the Particle Data Group (PDG) [H] average of the scale-invariant masses 
in the MS scheme [i.e. solutions to = /x]. 


mc{mc) = 1.275(25) GeV, 

(1) 

= 4.18(3) GeV, 

(2) 


is dominated by tuq extractions from low-energy data. These MS masses, as well as pole 
masses, have been used in the literature to estimate the theoretical precision achievable in 
precision Higgs calculations [iniiis]. 

However, looking into the future, such indirect engagement of low-energy observables in 
precision Higgs analysis might be ultimately unsatisfactory. A large amount of low-energy 
data has been highly processed to yield just two numbers, as in Eqs. Q and 0. It is not 
even clear whether these numbers accurately reflect our knowledge of mq, because the aver¬ 
aging involves mq extractions some of which are apparently correlated due to similar data 
and/or methods used. The error bars assigned to them contain experimental uncertainties 
from many different measurements, as well as theoretical uncertainties from calculating many 
different quantities. In addition, a self-described inflation of uncertainties by the PDG [16] 
is introduced to account for underestimated systematic errors in some mq extractions [T7] . 
Finally, Eqs. ([^ and (|^ do not retain possible correlations between as{mz) and the ex¬ 
tracted mq. They are thus treated as independent inputs in precision Higgs analysis, which 
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is strictly speaking not correct. 

As we strive for the highest-precision calculation possible in order to match percent (or 
even perhaps parts-per-mil) level of experimental precision achievable in the foreseeable 
future^, the rich information hidden in Eqs. ([^ and ([^ should be revealed, and the role of 
individual low-energy observables emphasized. Conceivably, a global £t would become 
more powerful in testing the SM when low-energy observables sensitive to mg as well as 
Higgs observables are incorporated. The scale-invariant masses mQ{mQ) would be then only 
inputs of the calculation. They are not considered as observables with experimental values 
and uncertainties, but are parameters to be tuned to minimize the function, where only 
true observables are included. 


In this paper we propose the idea of directly working with low-energy observables in pre¬ 
cision Higgs analysis. In addition to the global £t perspective mentioned above, low-energy 
observables can also play a role in identifying individual sources of theoretical uncertainties 
in precision Higgs calculations. This is conveniently done by eliminating mg (mg) from our 
input in favor of two low-energy observables, and recasting Higgs observables in terms of 
these and other input observables. For this procedure to be meaningful, the two observ¬ 
ables chosen should be representative of the large amount of low-energy data contributing to 
Eqs. 0 and ([^, in the sense that mg extracted from them alone should be precise enough. 
In the language of a global £t, the ideal choices would be two observables that dominate 
the low-energy observables contribution to y^. In this regard, a reasonable, though by no 
means exclusive, option would be to use the moments and of —)■ QQ inclusive 

cross section, dehned by 


/ ds 

—i?g(s), where Rg = 


a{e~^e~ 


—)■ QQX) 


(T(e+e“ —)■ /i+/i“) ’ 

with the precise dehnition of Rg from experimental data discussed in |22]- '^c(^c) and 
mb(mb) reported in the literature from analyzing these moments typically have O (10 MeV) 
uncertainties quoted [IIIE2H21]. For the Higgs observables we will focus on the partial widths 


^ Though precision measurements of Higgs observables, especially the partial widths into cc and bb discussed 
in this paper, are difficult at the LHC, such high precision is generally believed to be achievable at the 
International Linear Collider, the Future Circular Collider, and the Circular Electron Positron Collider. 
For recent analyses, see e.g. [181121) . We also note that for the bb channel, the importance of a higher 
theory precision is further emphasized by its relevance to the calculation of the total widths and all 
branching ratios of the Higgs boson. 
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Th^cc and T and assess the level of precision we can achieve in SM predictions for them. 
We will see that with direct contact made between these partial widths and the low-energy 
moments, the vagne notion of “nncertainties from tuq" is decomposed into concrete sonrces 
of nncertainties. In particnlar, parametric nncertainties from inpnt observables 
and as{mz) and pertnrbative nncertainties dne to missing higher-order corrections to the 
moments can be exposed separately. We note that while the parametric nncertainties are 
cnrrently expected to be at the percent level, and are in principle redncible with fntnre data 
and more carefnl experimental extraction of the moments, the pertnrbative nncertainties 
may represent a bigger challenge dne to lack of knowledge of the appropriate renormalization 
scales in the low-energy regime. It is therefore worthwhile to fnrther investigate theoretical 
as well as experimental aspects of the low-energy observables for the precision Higgs program 
to sncceed. 

II. INCORPORATING LOW-ENERGY OBSERVABLES INTO A GLOBAL 

PRECISION ANALYSIS 

The strongest tests of the SM rely on comparing its predictions across all accessible energy 
scales. By disentangling the information contained in the charm and bottom qnark masses 
in the context of precision Higgs analysis, we expose an interesting interplay between Higgs 
observables and low-energy observables. The sensitivity to mq that they share in common 
snggests the inclnsion of both in the precision program. 

An incomplete list of candidates for low-energy observables can be inferred from the 
mq extraction literatnre, and inclndes low [IZl I22ll24] and high [25H28] moments of Rq 
mentioned above, and their variants [29l [30], moments of lepton energy and hadron mass 
distribntions of semileptonic B decay pTHH3] . etc. We denote them collectively as {0)°'^}, 
with i rnnning from 1 to the nnmber of low-energy observables we wish to incorporate into 

^ It should be noted that we will treat asijnz) as both a calculational input and an observable with a central 
value and uncertainty. In principle one could treat as{mz) as merely a calculational parameter and let 
the observables that are highly sensitive to the as{mz) value be part of the global fit, analogous to what 
we have done with mq^mQ). However, as{mz) is one step further removed from direct determination of 
H bb, cc partial widths compared to mqimq), and so treating asijnz) as both an input parameter 
and (highly processed) observable is numerically justified. 
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the analysis. All these candidates should be carefully examined, and correlations among 
them should be understood, so that the best choices can be made for {0*°"}. 

In the high-energy regime, the observables include, for example, various partial widths, 
branching ratios, and production cross sections of the Higgs boson. Let us call them 
If not restricted to precision Higgs analysis, one may even include in the electroweak 

observables, such as the effective weak mixing angle, Z boson partial widths, and forward- 
backward asymmetries in e’''e“ annihilation at the Z pole. This will make the global analysis 
even more powerful, because the Higgs observables are sensitive to the same set of input 
observables as the electroweak observables: 

= {mz, Gf, a{mz), mt, as{mz), mn}- (4) 

Parenthetically we remark that the common practice of treating the top quark mass mt as 
an input observable is justihed for present purposes. A more careful treatment of mt, like 
what we do here with me and mb, may be needed in the future when precision measurements 
on the tt threshold are carried out at an collider. 

Additional calculational inputs, which are not necessarily of the observable type, include 
the charm and bottom quark masses {mg^mQ)} = {mdme), mb{mb)}■ There may be other 
input parameters, which we denote collectively by Examples are the r lepton mass, 

flavor angles, and nonperturbative parameters (e.g. gluon condensate) involved in some low- 
energy observables. 

Assuming the potentially complicated correlations among all the high- and low-energy 
observables will be understood in time, we may ultimately subject all the observables to a 


global fit, by minimizing the y 

^ function with respect to the inputs: 


Calculation inputs: 


(5) 

Fit observables: 

{0.} = {0”} u {ot**-} u {o;“}, 

(6) 

To minimize: 

x" = Elof'dii}) - or‘]y-‘[of({it}) - Of- . 

(7) 


ij 


Here “th” and “expt” denote theoretical and experimental values, respectively, and V is the 
covariance matrix containing uncertainties and correlations among observables. The calcu¬ 
lational inputs could just as well be chosen to be a minimal set of Lagrangian parameters; 
however, it is most convenient for our purposes to choose a combination of observables and 
Lagrangian parameters as the minimal set of calculational inputs. 
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Compared with the conventional approach where low-energy data contribute indirectly 
via the averaged {mg(mg)}, our proposal of directly working with low-energy observables 
allows appropriate treatment of all the correlations and uncertainties. In particular, there 
is no averaging over correlated mg extractions, and the calculational inputs {'mQ{mQ)} and 
o^si'^z) are no longer correlated. Challenging as it is, such a global analysis is worth further 
investigation. As a long-term goal for the precision program, it will test our understanding 
of elementary particle physics at an unprecedented level. 

As a hnal remark in this section, the techniques described above are to be employed in a 
rigorous test of the SM. The resulting statistical test from the analysis is for determining 
the likelihood of the compatibility of the data with the SM hypothesis. It is straightforward 
to apply these techniques to a slightly different model, which we call the kSM, dehned to 
be exactly the SM theory except that each coupling of the Higgs boson to SM states has a 
free parameter k* in front that is varied to £t the data (see e.g. [miallEH]). In that case, 
the analysis must include these k* as extra input variables and the resulting £t tests the 
compatibility of the kSM theory with the data and, if compatible, gives confidence intervals 
for the Ki values. Just as with the SM, at the next level of precision analysis of the kSM it 
is important to address the role of low-energy observables that we study in this paper. 


III. RECASTING HIGGS OBSERVABLES IN TERMS OF LOW-ENERGY 

OBSERVABLES 


In order to investigate sources of theoretical uncertainties in calculating the Higgs ob¬ 
servables, it is helpful to recast them in terms of a set of input observables without invoking 
a global £t. In the simplest case, suppose all the observables under consideration are in¬ 
sensitive to We choose two low-energy observables 0}°’", By inverting the 

functions 




Qlow ^ Q\ow 


{0“},{mg(mg)} 


( 8 ) 


we express the quark masses in terms of 0^°", 02 °*: 


rricimc) = {0“}, Of", 02°* , = m^(rrn,) {Of}, Of", 0'°" 


(9) 
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{mQ(mQ)} can then be eliminated from the calculation of the Higgs observables: 




^low 


( 10 ) 


and we have achieved the goal of recasting Higgs observables in terms of low-energy input 


observables Oj”", 0 ^ 2 '^. From Eq. (10) it is clear that the precision in the SM prediction for 


the Higgs observables will benefit from improved knowledge of mg, which ultimately comes 
from better measurements of the low-energy observables. 

Our choices for the low-energy input observables, 


= Ml Ml 


( 11 ) 


require only a slight generalization of the simple formalism above. We will take into account 
an additional input, the gluon condensate, as in the case of Ml, but its contribution 

allows for a simplihed treatment. In fact, the simplicity of the analysis is our main motivation 
for choosing these moments as inputs rather than other low-energy observables which lead 
to similar level of precision in the extracted mg. For example, if we were to use semileptonic 
B meson decay observables (see e.g. [niH33]), more input parameters in will show 

up, including flavor angles and four nonperturbative parameters. Also, the low moments 
{M^ with 77, < 4) chosen here are computationally more straightforward than the high 
moments (n > 10; see e.g. [251 - 128] ). The former can be calculated conveniently in the 
relativistic theory, while a nonrelativistic effective theory treatment is needed for the latter. 
In addition, since the calculation involves MS quark masses, there is no need for introducing 
other mass schemes. Potentially large uncertainties associated with mass scheme conversion 
(e.g. from pole or kinetic masses to MS masses), which is needed for some other methods, 
can thus be avoided. We also note that the approach of extracting mg from the low moments 
was recently recast by the lattice QCD community |36H38] . and future development in this 
direction may shed light on the precision Higgs program [3^ . 

To calculate one applies quark-hadron duality mi to relate the moments Ad^ to 
vector current correlators. 





where 




( 12 ) 


(13) 









with being the electromagnetic current of Q. Hq can be calculated as an operator product 


expansion: 


■W? = i: Irf + X?™. (14) 


where Qq is the electric charge of quark Q. As one can see, the values of these moments 
depend on the quark masses, a fact that QCD sum rules practitioners use to extract quark 


masses (for reviews see [HI Il2])- The two terms in Eq. (14) come from perturbation the¬ 
ory and nonperturbative condensates, respectively. The perturbative part is known up to 
O {a^) [13], while the gluon condensate contribution, which dominates has been 

calculated to next-to-leading order M- Note that the coefficients are functions of n/, 
the number of active quark ffavors. The common choices are n/ = 4 for Q = c and n/ = 5 
for Q = b. These are also the numbers of active quark flavors assumed for as{fi) and mQ{fi) 


in Eq. (14). as{mz) is defined for uj = 5, and should be matched to the 4-flavor effective 
coupling at the bottom quark threshold before being used in Eq. (14) for Ai^- 
culations the matching is done assuming 4.2 GeV for both the threshold scale and 
but all the results are found to be insensitive to the details of threshold matching. 

mQ{fi) are usually extracted by comparing the theoretical calculation with experimental 
data for (see [IZ1122I for technical details). Normally the lowest moment Ai^ is taken 
for the charm quark so as to suppress the nonperturbative contribution to the subpercent 
level Ha Ea US). For the bottom quark the gluon condensate can be safely neglected at 
the present level of precision [22], and the second moment Ad 2 is preferred due to large 
experimental uncertainty in Ad^ We also neglect O {rril/ml) terms in Ad^, not explicitly 


written out in Eq. (14), which constitute a tiny contribution 

It is pointed out in ini that the scales at which tuq and as are renormalized should be 


considered independently to avoid bias in the uncertainty estimate. Eq. (14) then should be 
generalized to 


kjQ _ (Qq/( 2/3)) ^(a,fe)/ X b 

{2mQ{fim)) i,a,b ^ ^ 


In' 




K 


K 


+ Ad«’"P. (15) 


The coefficients in this equation can be readily derived from via renormalization 

group (RG) equations, and numerical results foru/ = 4 can be found in [T7]. Due to unknown 
O (af) terms, the calculated Ad(^ exhibit dependence on both and pa. Scale dependence 
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is a general feature of finite-order perturbative calculations, and should be considered with 
care in estimating theoretical uncertainties. We have more to say on this below. 


With mQ(/im), C(s{fJ-a) related to mQ^mg)^ as{mz) via RG equations, Eq. (15) matches 
the general form of Eq. ([^, with agimz) being the only relevant element in {O™}. There 
are additional inputs /^a and So in our case, Eq. ([^ is modihed as: 

= M\ y ( 16 ) 

Ml = Ml[as{mz),mb{mb),fil^,fil], (17) 

where we have neglected M 2 ^^- As mentioned above, the nonperturbative contribution has 
been claimed to be negligible for the bottom quark. We have checked this in the case of 
Ml, where the contribution from MI’'^^ is below 0.1%, which should be compared to the 
experimental uncertainty of Ml of about 1%. Treating Mi^'^ and mc{mc) as independent 
inputs, which we will justify later, and focusing on the Higgs boson partial widths to cc and 


bb as examples of we have, in place of Eqs. ([^ and (10), 

mcimc) = mcimc)[as{mz), Ml, 
mb{mb) = mb{mb) as{mz). Ml, , 

hiJ—^cc } y y 

^ H^bb = ^ H^bb 


(18) 

(19) 

( 20 ) 


= T 


H^bb 


( 21 ) 


{OT}ymb{mb),HH 

where collectively denote other renormalization scales involved in the calculation of 

the partial widths. These are nevertheless not the only scale dependences for the partial 
widths in such an analysis. The residual scale dependences of the low-energy observables 
are seen to propagate into the extracted quark masses, and constitute part of the uncer¬ 
tainties in rngirriQ). These uncertainties eventually propagate into the calculations of Higgs 


observables, and are reflected in the /Xm, dependences in Eqs. (20) and (21). Note also 


that in the second equalities in Eqs. (20) and (21), the as{mz) dependence in the partial 


widths has been changed to account for the correlation with rngimo) reflected in Eqs. (18) 


and (19). 


Eqs. (20) and (21) represent the hnal results of the exercise of recasting Higgs observables 
in terms of low-energy observables, with the information contained in mg (mg) fully resolved. 
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They will be used in the next section to investigate the theoretical uncertainties in these 
partial widths. 

To close this section we remark on the treatment of The known terms read 


= 




(2mr'y L 


-16.042 - 168.07 




+ 


TT 




( 22 ) 


where is the gluon condensate. The commonly used value in the context of charm 

quark mass extraction is derived from r decay data 


a.. 


TT 


—G" ) = 0.006 ± 0.012 GeVh 


(23) 


In addition to the imprecise knowledge of ^^G^y, we note two other sources of uncertainties 
in First, it is argued in [T^ SS] that should be expressed in terms of the pole 

mass rather than the MS mass in order to have a stable a* expansion. We agree with 
this argument, but note that the use of the pole mass may introduce further ambiguities. 
For example, if one tries to calculate the pole mass from the MS mass, the result will be 
very sensitive to the loop order. Second, considerable uncertainty is introduced by the 
fi dependence of the bracket in Eq. (22), since the O (a^) terms are not known. This 


renormalization scale is not necessarily related to or in fhe perturbation theory 


contributions [the hrst term in Eq. (15)]. All these uncertainties and ambiguities will dilute 
any conceivable correlation between and mc{mc), justifying our treatment of them as 

independent inputs. In our analysis the following value for will be assumed: 


^c,np ^ _0.000i;[]:[][]?^ GeVl 


(24) 


The central value corresponds to = 0.006 GeV'^, = 1.7 GeV and /r = 3 GeV 


in Eq. (22). The errors are very conservatively estimated by taking the extreme values 


^poie _ GeV, jjL = 1 GeV, and varying in the range in Eq. (23). Even with the 

extreme values considered, is still a subpercent-level contribution to ~ 0.2 GeV“^. 


IV. THEORETICAL UNCERTAINTIES OF HIGGS PARTIAL WIDTHS 


It is clear from Eqs. (20) and (21) that there are two types of uncertainties in the calcu¬ 


lation of the Higgs partial widths. Parametric uncertainty results from imprecise knowledge 
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of the input parameters, including the input observables (A^i, those in {0™}) 

and the nonperturbative parameter The experimental values and errors of the input 

observables are: 


Ml 

= 0.2121(20)(30) GeV-2 [I7|, 

(25) 

M\ 

= 2.819(27) X IQ-^GeV^ |15], 

(26) 

as{mz) 

= 0.1185(6) [m. 

(27) 

rriH 

= 125.7(4) GeV [H], 

(28) 

mt 

= 173.21(51)(71) GeV [H], 

(29) 

mz 

= 91.1876(21) GeV [H], 

(30) 

a{mz) 

= 1/127.940(14) HU, 

(31) 

Gf 

= 1.1663787(6) x 10"® GeV^ [Hj. 

(32) 


For A4f and mt the two experimental uncertainties are statistical and systematic, respec¬ 
tively. There is an additional systematic uncertainty in AI 2 associated with the prescriptions 
used in extracting moments from data. This is discussed in [33], and we adopt “Option A” 
in that paper because among the three options considered there it appears to yield the most 
consistent results for mg (mg) across different moments. 


Perturbative uncertainty, on the other hand, is associated with unknown higher-order 
terms in perturbation theory calculations, and leads to residual dependence of calculated 
observables on the renormalization scales. When the partial widths are recast in terms of 


A4f and as in Eqs. (20) and (21), multiple scales enter, fin comes from the calculation 
of the Higgs boson decay. The associated perturbative uncertainty has been studied in the 
literature; see e.g. [10] where it is found to be small compared with parametric uncertainty. 
Here we focus on /i^, /io, which originate from the calculation of the low-energy observables 


M 2 [see Eqs. (15 17)]. Their contribution to the total theoretical uncertainty will be 


singled out below by setting all input parameters to their central values in Eqs. (24 32), and 
setting jUff = mu- 


We study the perturbative uncertainty from /i^, /Uq, in two steps. First, mg^firn) are 


calculated by iteratively solving Eq. (15) following the procedure explained in |IZ|. from 


which mQ{mQ) are derived. We use the RunDec package |37| for RG running and threshold 
matching to the highest loop order implemented in the package. Second, the partial widths 
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FIG. 1: Contours of mc{mc) (top-left), mb{mb) (bottom-left) in GeV, and Th^cc (top-right), T 
(bottom-right) in MeV in the plane. These plots demonstrate Eqs. ( TspI ) with all other 

inputs hxed. The unlabeled contours represent decreasing values toward the top-left corner in steps 
of 0.01 GeV, 0.005 GeV, 0.002 MeV, 0.005 MeV, respectively. 


Th^cc, ^H^bb calculated using the expansion formulas in |T0]. The results of both steps 

are shown in Fig. [^as contour plots in the plane^. They correspond to Eqs. (18 21) 


^ The numerical difference between our rricijnc) contour plot and Fig. 6(c) in |17j is due to the input 
and as{mz) used, and to a lesser extent the treatment of At 
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with other inputs fixed. These plots illustrate the propagation of /im, dependence from 
low-energy moments calculations to Higgs partial widths. 


To estimate the perturbative uncertainty, a common practice is to identify a characteristic 
scale of the process of interest, and vary the renormalization scale within a factor of two 
around that scale. For example, has been varied from mH/‘i to 2mH in [lU]- However, 
this method is not directly applicable to and /Iq,, since receive contributions from 
all energy scales as evident in Eq. (|^. One might guess from qualitative features of Rq{s) 
that the characteristic scale should be O (2mQ), the masses of quarkonium resonances. But 
due to the relatively large value of a* in the low-energy regime, the exact number, and hence 
the range in which we choose to vary /i^, /^a can greatly affect the result of our uncertainty 
estimates. This is already clear from Fig. where Th^cc and are seen to exhibit 

rapid variation in the low-/rm regime. 


Lacking an optimal method to estimate the perturbative uncertainty, we refrain from giv¬ 
ing exact numbers, but instead aim to illustrate the ambiguity in the estimate of perturbative 
uncertainty by varying pm and Ha independently within an adjustable range [pmin, We 

will focus on the uncertainties in the partial widths, and remark that they are related to the 
uncertainties in mQ^mq) by [10] 


AF 


H^cc 


^rncjmc) ^ ^ H^bb _ 


Amfc(mfe) 


X 0.56%. 


(33) 


^H^cc 10 MeV T 10 MeV 

The perturbative uncertainty, defined as half the difference between the maximum and min¬ 
imum values of F/^^cc, ^H^bi^ depends on pmin and Pmax- We present the results in Fig. 
in terms of “percent relative uncertainties,” defined to be lOOAF/F. The red solid curves 
show the estimated perturbative uncertainties as functions of /imin, with /rj^ax (/^max) fixed 
at 4 (15) GeV. Alternative choices for /ijnax (/^max); 3 and 5 (13 and 17) GeV, give rise to 
the red dashed curves. These can be compared with the dominant parametric uncertainties 
shown by the other curves in Fig. (see figure caption for details). The popular choices in 
the literature (/x^ia,/x^ax) = (2,4) GeV and (/imin,FLx) = (5,15) GeV yield perturbative 
uncertainties of 1.2% and 0.33% for Th^cc and F n^bb^ respectively, comparable with para¬ 
metric uncertainties. However, the perturbative uncertainties increase rapidly and dominate 
the total theoretical uncertainties if lower renormalization scales are considered. The result 
of the theoretical uncertainty estimate is then strongly dependent on the artificial choice of 
fimin- This poses a serious ambiguity in precision analysis, and calls for more enlightened 
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FIG. 2: Percent relative uncertainties in Vh^cc (left) and T (right) as functions of /imin 
from various sources: perturbative uncertainty with //(nax = ^ GeV, GeV (red solid) 

or alternatively /u(^ax = 3,5 GeV, /x^ax = 13,17 GeV (red dashed), parametric uncertainties 
from Jv[\ or (orange), as{mz) (cyan solid), (blue, for Fh^cc only) and mn (purple). 

The parametric uncertainty from as{rnz) incorrectly calculated assuming no correlation with mg 
(cyan dotted) is also shown for comparison. The parametric uncertainties are defined as shifts of 
the central values of Th^cc and T for ^min < < ^max caused by varying the input 

parameters within the errors quoted in Eqs. ( 24||^ ), with = 4 GeV, /x^ax = 15 GeV (the 
kinks are due to the maximum or minimum shifting to a different region in the plane), and 

are found to be insensitive to ^max- 


prescriptions for the uncertainty estimate. We note two possible directions in this regard. 

The hrst direction was suggested very recently in |18] in the context of mg extraction. 
There it is argued that the large perturbative uncertainty from completely uncorrelated 
variation of /x^ and /Xq, is probably an overestimate. To get the perturbative uncertainty 
under control, a “convergence test” is performed to identify regions in the /Xm-ho plane 
where the perturbative series converges too slowly (characterized by a large convergence 
parameter). These regions are then discarded in the uncertainty estimate. Following the 
approach outlined in |18] , we find that the discarded regions correspond to the upper-left and 
bottom-right corners in each plot in Fig. where mg{mg) and the partial widths exhibit 
rapid variation. The final result in |1H] is a reduced perturbative uncertainty: 14 MeV and 
10 MeV for mcirric) and mh{mb), respectively, corresponding to 2.9% and 0.56% relative 
uncertainties in Th^ce and Frespectively. 
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The convergence test is a well-motivated idea, reflecting the intuition that a proper scale 
choice should not lead to very slow convergence. However, further study is necessary to 
examine various details of the approach. For instance, one may consider loosening the con¬ 
straints mcirric) < — 4 GeV, mhijnb) < lAi — 15 GeV imposed in [IH]. In 

particular, pm, slightly lower than rngimo) should be allowed as long as one retains 
4-flavor (5-flavor) effective strong coupling for the charm (bottom) quark. Also, the conver¬ 
gence criterion may be refined. The definition of the convergence parameter in |18] assumes 
an approximate geometric series behavior of the series, but we And the latter falls off 
more slowly than a geometric series in most cases. Furthermore, it remains to seek a less 
arbitrary prescription for the fraction of (/Xm, fia) to be discarded, and to investigate whether 
the convergence parameter is a good indicator of the size of higher-order corrections. In any 
case, to be conservative the reduced perturbative uncertainties mentioned above should be 
interpreted with caution before the approach is developed further. 

As an alternative direction, one may consider the possibility of finding an optimal scale 
via a defensible scale-setting procedure, such as the one advocated by Brodsky-Lepage- 
Mackenzie (BLM) |19]. The BLM scale for an observable is obtained by absorbing the 
rif terms in the perturbation series, which come from the QGD beta function, into the 
running coupling This is arguably the physical scale of the process, with higher-order 
corrections associated with RG running appropriately resummed. We also note that the 
BLM procedure extended to all orders based on the principle of maximum conformality [50] 
has been demonstrated to be self-consistent [51]. In the case of Ad®, however, there are two 
renormalized parameters and mg, and naive application of the BLM procedure might be 
problematic. This is because even when the nj terms are absorbed into running ag and/or 
mg, the leading-order mass renormalization, which is independent of nj, may lead to large 
loop corrections which are difficult to identify. Indeed, we find that naive application of BLM, 
namely absorbing the Ufol terms, sets scales for pm and /Iq, which are strongly disfavored 
by the convergence test. In light of the importance of a more precise mg determination, it 
might be worthwhile to investigate the nontrivial possibility of generalizing the BLM method 
and its extensions [501150] to include running quark masses. 

The parametric uncertainties, on the other hand, are seen from Fig. [^ to be dominated by 
experimental measurement uncertainties of and Ai\ (orange). Reduction of these will 
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rely on more precise measurements of Rq{s) and more careful treatment of experimental 
data. At present the major problem is the lack of data above y/s = 11.2 GeV, resulting in 
large uncertainties in the bottom quark moments [H]. Also, the quarkonium resonances are 
currently treated in the narrow width approximation, the quality of which should be exam¬ 
ined in light of higher precision requirements in the future. as{mz) (cyan solid) constitutes a 
subdominant source of parametric uncertainties. Its contribution is seen to be smaller than 
the incorrect estimate assuming no correlation between as{mz) and tuq (cyan dashed), due 
to partial cancelation between direct as{mz) dependence and indirect dependence through 


rriQ. With our conservative estimate (i.e. erring on the large side) in Eq. (24), leads to 

an uncertainty in Th^cc (blue) at a similar level as as{mz)- This may represent a challenge 
in the future, and calls for further investigation of the gluon condensate contribution. The 
uncertainty due to mu (purple) is less important, while other input observables listed at the 
beginning of this section have a negligible effect on the parametric uncertainty. 


V. CONCLUSIONS 

For the precision Higgs program to succeed in the future, additional effort is required to 
improve the precision of SM calculations in order to match the proposed experimental ac¬ 
curacy. A better understanding of theoretical uncertainties is critical. Toward this aim, we 
emphasize the role of low-energy observables, and further propose the idea of a global anal¬ 
ysis incorporating relevant observables across all energy regimes. Rather than contributing 
indirectly via the charm and bottom quark masses, low-energy observables explicitly par¬ 
ticipate in such a precision analysis. Future studies in this direction should examine all 
candidates of low-energy observables, and determine an efficient set of observables for the 
global fit. 

In the context of precision Higgs calculations, we focused on the Higgs boson partial 
widths to charm and bottom quarks, and investigated the theoretical uncertainties in these 
observables. By eliminating charm and bottom quark masses in favor of low-energy observ¬ 
ables }A\ and we recast the partial widths in terms of these and other input observables. 
Much information originally hidden in uncertainties in the highly processed quark masses 
becomes transparent. Experimental uncertainties in the low-energy observables are directly 
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propagated into the Higgs partial widths, and the uncertainty due to asijnz) is treated 
properly. Perturbative uncertainties are difficult to assess due to the ambiguity in the choice 
of renormalization scales in the low-energy regime, and can dominate the total theoreti¬ 
cal uncertainty of the Higgs partial widths if lower values of the renormalization scales are 
considered than is usually the case in the literature. 

Such analysis points to future directions in the precision program. For the partial widths 
considered here, we note that while future experimental progress could potentially reduce 
parametric uncertainties signihcantly, our ability to make precise predictions on the Higgs 
partial widths will not improve unless better understanding of the perturbative uncertainty 
is achieved. As for and M.\ studied here, this might require the calculation of O (a^) 
corrections to nQ(g^) (in the low-g^ limit) and/or more enlightened scale setting. Though the 
actual situation may be better in a global £t where Ai\ and Ai\ the only low-energy 

observables involved, it remains crucial to carefully investigate whether the scale-setting 
problem is also present for other low-energy observables sensitive to mq. If the perturbative 
uncertainty gets under control, the precision program, where both low-energy observables 
and Higgs observables play an important role, will be promising in studying properties of 
the Higgs boson, and even more generally testing the SM across a wide range of energy 
scales and probing new physics ideas. 
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